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In this paper we study the H + 7 (H — and T) production in coherent hadron - hadron 

interactions at LHC energies. Considering the ultrarelativistic protons as a source of photons, we 
estimate the 7+p — > H + 7 + X cross section using the non-relativistic QCD (NRQCD) factorization 
formalism and considering different sets of values for the matrix elements. Our results for the total 
p + p — >p + H + j + X cross sections and rapidity distributions at ^fs = 7 and 14 TeV demonstrate 
that the experimental analysis of the J/$ + 7 production at LHC is feasible. 
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I. INTRODUCTION 



In the last years, the analysis of coherent hadron-hadron collisions becomes an alternative way to study the theory 
of strong interactions - the Quantum Chromodynamics (QCD) - in the regime of high energies (For reviews see 
Ref. i]). The basic idea in coherent hadron collisions is that the total cross section for a given process can be 
factorized in terms of the equivalent flux of photons in the hadron projectile and the photon-photon or photon- 
hadron production cross section. The main advantage of using colliding hadrons and nuclear beams for studying 
photon induced interactions is the high equivalent photon energies and luminosities, that can be achieved at existing 
and future accelerators. Consequently, studies of jp interactions at LHC provides valuable information on the QCD 
dynamics at high energies. The photon-hadron interactions can be divided into exclusive and inclusive reactions. In 
the first case, a certain particle is produced, while the target remains in the ground state (or is only internally excited). 
On the other hand, in inclusive interactions the particle produced is accompanied by one or more particles from the 
breakup of the target. The typical examples of these processes are the exclusive vector meson production, described by 
the process jp — > Hp (H = J/ 1 ^, T), and the inclusive heavy quark production [jp —> XY (X — cc, bb)], respectively. 
The results of these studies demonstrate that their detection is feasible at the LHC (For recent discussions see, e.g. 
Refs. @-(3])- It motivates the analysis of the production of other final states in coherent hadron - hadron interactions. 

In this paper we study, for the first time, the inclusive quarkonium + photon photoproduction in pp collisions at 
LHC energies. The total cross section for the process p + p — >p<E>H + ^ + X (H = J/^> or T) is given by 

a(pp^p®H + <y + X) = 2 J du> a lp ^H+ 1+ x{W^ p ) , (1) 

where <g> represents a rapidity gap in the final state, u is the photon energy, is the equivalent photon flux, 
W^ p = 2 wVSnn and V^nn is the c.m.s energy of the hadron-hadron system. The photon spectrum of a relativistic 
proton is given by 0, 
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with the notation = 1 + [ (0.71 GeV 2 )/Q^ in ], Qf nin = uj 2 /[-fl (1 - 2u;/ v / 3W)] » (w/ 7L ) 2 and 7L is the Lorentz 
boost of a single beam. Distinctly from the exclusive vector meson photoproduction discussed in, e.g., Ref. [3j, which 
is characterized by two rapidity gaps in the final state, in our case it is characterized by one rapidity gap associated to 
the photon emitted by one of the protons and the remnants of the other proton. The main input in our calculations 
is the 7 + p — > H + 7 + X cross section, which we estimate using the non-relativistic QCD (NRQCD) factorization 
formalism Q. Recent results demonstrate that this formalism is able to describe quite well the RHIC and HERA 
data for the hadroproduction and photoproduction of charmonium in terms of an universal set of matrix elements 
and that the inclusion of the color octet processes is indispensable to describe the photoproduction data [3] ■ 

This paper is organized as follows. In next section we present a brief review about the quarkonium+7 photopro- 
duction in the NRQCD formalism. In Section IlIII we present our predictions for the rapidity distributions and total 
cross sections for J/'I' + 7 and T + 7 production at LHC energies. Finally, in Section llVl we summarize our main 
conclusions. 
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FIG. 1: Energy dependence of the total J/^ + 7 photoproduction cross section considering two different sets of matrix elements 
(see text). 



II. THE QUARKONIA+7 PHOTOPRODUCTION 



In the NRQCD formalism the cross section for the production of a heavy quarkonium state H factorizes as a(ab — > 
H + X) — J2„ ^{o-b ^ QQ[n] + X)(O h [n]), where the coefficients a(ab —> QQ[n] + X) are perturbatively calculated 
short distance cross sections for the production of the heavy quark pair QQ in an intermediate Fock state n, which 
does not have to be color neutral. The (0 H [n\) are nonperturbative long distance matrix elements, which describe 
the transition of the intermediate QQ in the physical state H via soft gluon radiation. Currently, these elements have 
to be extracted in a global fit to quarkonium data as performed, for instance, in Ref. Q. It is important to emphasize 
that the underlying mechanics governing heavy quarkonium production is still subject of intense debate (For a recent 
review see, e.g., Ref. @). 

In the specific case of the H + 7 photoproduction, the total cross section can be expressed as follows [lol ] 



cr(7+p^ H + j + X) 
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where z = (j>h -p)/(p-y-p), with p and p 7 being the four momentum of the quarkonium, proton and photon, 
respectively. In the proton rest frame, z can be interpreted as the fraction of the photon energy carried away by 
the quarkonium. Moreover, p± is the magnitude of the quarkonium three-momentum normal to the beam axis. The 
partonic differential cross section da/dt is given by [l(| 



do , . 
— (7 + .9-> H + 7) = 



64tt 



eQa 2 a s m Q 



s 2 s\ + t 2 t\ + u 2 u\ 



s\t\u\ 



(4) 



where eQ and tuq are, respectively, the charge and mass of heavy quark constituent of the quarkonium. The Man- 
delstam variables can be expressed in terms of z and p± as follows: 



pl + (2m Q ) 2 (l-z) 
z(l-z) 
p\ + (2m Q ) 2 (l - z) 



l-z 



(5) 
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FIG. 2: Rapidity distribution for the +7 photoproduction in pp collisions at y/s = 7 TeV (left panel) and 14 TeV (right 
panel) considering two different sets of values for the matrix elements. 



Moreover, s\ — s — 4toq, t\ = t — 4wiq, u\ = u — 4m q and the Bjorken variable x can be expressed by 

= pj + (2m Q ) 2 (l-z) 
W* p z(l - z) ' 



(6) 



where W lp is the photon - proton center-of-mass energy. In our calculations we consider the gluon distribution as 
given by the CTEQ6 parametrization [Tlj . 



III. RESULTS 



In what follows we present our results for the H + 7 photoproduction in coherent hadron - hadron collisions at LHC 
energies. One of the main uncertainties in our predictions is associated to the NRQCD matrix elements. For the J/ty 
case, we consider the more recent global fit performed in Ref. [1] (denoted BK hereafter), and its uncertainties, which 
obtain that (0^ N '{ 3 Si)) = 2.24 ± 0.59 x KT 3 GeV 3 . For comparison, we also use the value previously considered 

in Ref. [T^ (denoted MEHEN hereafter): (0^( 3 Si)) = 6.6 x 1CT 3 GeV 3 . Although this choice is outdated, we 
keep it in order to estimate the dependence of the matrix elements in our results for the J/'F + 7 production in 
coherent interactions. As show in Fig. [TJthese two choices imply very distinct results for the J/\f/ + 7 photoproduction 
cross section, with the latter being an upper bound for the total cross section. The cross section increases with the 
energy, which is directly associated to the x-behaviour of the gluon distribution. It is important to emphasize that 
while studies of photoproduction at HERA were limited to photon-proton center of mass energies of about 200 GeV, 
photon-hadron interactions at LHC can reach one order of magnitude higher on energy. For instance, if we consider 
pp collisions at LHC, the Lorentz factor is jl — 7455, giving the maximum c.m.s. 7p energy W JP s» 8390 GeV. 
Consequently, the analysis can also be useful to constrain the proton gluon distribution. (For related studies see Ref. 
[l2l]). We postpone a more detailed study of this topic for a forthcoming publication. 

For the T case, we also use two different sets of values for the matrix elements. The first one taken from Ref. [l3[ 
(denoted BSV hereafter), which obtain (Oj( 3 Si)) = 5.3 x 10~ 3 ± 0.5 GeV 3 from fits to CDF data for bottomonium 
production. The second one we taken from Ref. [14[ (denoted BFL hereafter) which obtain (Oj( 3 Si)) = 0.02 GeV 3 . 

Lets now calculate the rapidity distribution and total cross sections for H + 7 production in coherent proton - 
proton collisions at -y/s = 7 and 14 TeV. The distribution on rapidity Y of the produced final state can be directly 
computed from Eq. (|TJ) , by using its relation with the photon energy oj, i.e. Y oc In {oj/mu). Explicitly, the rapidity 
distribution is written down as, 

da\ P + p^p®H + 1 + X)\ dN l/hl {uj) dN l/h2 (uj) 

dY = ^ du) a lh2^H+ 1 +X (U) + UJ — (T-fhi-W+j+X { UJ ) i I'J 

where we taken into account that the two protons can be the source of the photons ( hi = h 2 = p) and X is a hadronic 
final state resulting of the fragmentation of the proton target. The rapidity gap ® is associated to the proton which 
emits the photon and remains intact. In Fig. [5]we present our results for the J/^ + 7 production considering the two 
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TABLE I: The total cross section for the H + 7 photoproduction in coherent hadron - hadrons collisions at LHC energies. 




FIG. 3: Rapidity distribution for the T ± 7 photoproduction in pp collisions at y/s = 14 TeV considering two different sets of 
values for the matrix elements. 



different sets of values for the matrix elements and two different center-of-mass energies. As the photoproduction cross 
section is proportional to the cross section, our predictions for the rapidity distribution increases with the energy. 
Moreover, as expected from Fig. [TJ the predictions obtained using the MEHEN matrix element [10( implies a value 
for rapidity distribution at Y — which is a factor w 3 larger than that using the BK one. This large difference is also 
observed in the predictions for the total cross section (See Table J]). It demonstrates that the study of this observable 
can be useful to constrain the matrix element. Assuming the design luminosity C = 10 7 mb _1 s _1 the corresponding 
event rates will be sw 1 x 10 5 (2 x 10 5 ) events/years at yfs = 7 (14) TeV for the BK matrix elements. 

In Fig. [3]and Table|I]we present our predictions for the T±7 production considering the two different sets of values 
for the matrix elements and yfs = 14 TeV. As observed in the J/vP case, the predictions are strongly dependent on 
the matrix elements used in the calculations. However, in the T case the cross sections are a factor 10 3 smaller than 
those obtained for J '/\V +7 production, which implies that the event rates will bewlx 10 2 events/years, making the 
experimental analysis of this final state a hard task. 

Lets discuss the experimental separation of the H + 7 photoproduction and compare our predictions with those 
obtained for the inclusive production or considering diffractive (Pomeron) interactions. In comparison to the inclusive 
hadroproduction (See e.g. [15]), which is characterized by the process p ± p — » X + H + j + Y 7 with both proton 
producing hadronic final states, the photoproduction cross sections are a factor fs 10 3 smaller. However, as in 
photoproduction we have a rapidity gap in the final state, the separation of the signal from hadronic background 
would be relatively clear since the event multiplicity for photoproduction interactions is lower, which implies that it 
may be used as a separation factor between these processes. Another mechanism which also is characterized by one 
rapidity gap in the final state, is the single diffractive H + 7 production. This process was analysed, e.g.. in Refs. 
[Ha . Il7| . considering that the Pomeron has a partonic structure. In comparison with those results, our predictions 
are a factor 8 smaller. However, it is expected that emerging hadrons from Pomeron processes have a much larger 
transverse momentum than those resulting from photoproduction processes. Consequently, in principle it is possible 
to introduce a selection criteria to separate these two processes. Moreover, it is important to emphasize that Pomeron 
predictions are strongly dependent on the value used for the gap survival factor, while our results should not be 
modified by soft absorption corrections. 
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IV. CONCLUSIONS 

In this paper we have computed for the first time the cross sections for photoproduction of quarkonium+7 in 
coherent pp collisions at LHC energies using the NRQCD formalism and considering two different sets of values for 
the matrix elements. Such processes are interesting since the final state is characterized by a low multiplicity and one 
rapidity gap. Moreover, the produced large px quarkonia are relatively easy to detect through their leptonic decay 
modes and their transverse momenta are balanced by the associated high energy photon. Our results demonstrate that 
the rapidity distributions and total cross sections are strongly dependent on the magnitude of the matrix elements. 
Moreover, we predict sizeable for the J/* + 7 cross section, which makes the experimental analyses of this process 
feasible at LHC. 
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